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We theoretically investigate the collective modes of imbalanced two component one-dimensional 
Fermi gases with attractive interactions. This is done for trapped and untrapped systems both at 
zero and non-zero temperature, using self-consistent mean-field theory and the random phase ap- 
proximation. We discuss how the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state can be detected 
and the periodicity of the associated density modulations determined from its collective mode spec- 
trum. We also investigate the accuracy of the single mode approximation for low-lying collective 
excitations in a trap by comparing frequencies obtained via sum rules with frequencies obtained 
from direct collective mode calculations. It is found that, for collective excitations where the atomic 
clouds of the two spin species oscillate largely in phase, the single mode approximation holds well 
for a large parameter regime. Finally we investigate the collective mode spectrum obtained by 
parametric modulation of the coupling constant. 
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I. INTRODUCTION 

With the advancements in the field of ultracold atomic 
gases new phases of matter are becoming experimentally 
accessible. Superfluidity [l|, [I] and density imbalance 
01 in two component Fermi gases have been demon- 
strated experimentally. Together these offer the ex- 
citing perspective of achieving the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase 00], an unconventional su- 
pcrfluid state in imbalanced Fermion systems. It has 
been under theoretical investigation since the 1960s but 
has so far eluded unambiguous experimental detection. 

In solid state systems the detection of the FFLO phase 
is unfavourable since density imbalances between spin up 
and down electrons need to be imposed via Zeeman split- 
ting, requiring strong magnetic fields. This in turn de- 
stroys superconductivity via diamagnctic effects Q. Ul- 
tracold atomic gases arc favourable for its detection since 
it is possible to set the density imbalance of the system 
simply by adjusting the number of spin up and down 
particles which are loaded into the atomic trap 0, 3 ■ 

Theory does predict the occurrence of the FFLO phase 
in three-dimensional (3D) ultracold Fermi gases though 
the region of parameter space it occupies is thought to be 
relatively small (for a review see [1]). On the other hand 
the FFLO phase is expected to be stabilised in quasi one- 
dimensional (IE)) geometries which can be imposed using 
optical lattices fllMil- 

In a recent development one-dimensional imbalanced 
Fermi gas systems have been created in experiment [isj . 
These recent experimental results indicate that the den- 
sity distribution is of the form expected for the partially- 
polarised state associated with the FFLO phase. How- 
ever, a direct demonstration of the ordering in the FFLO 
phase remains to be shown. An important question which 
now arises, is how best to detect this characteristic or- 
dering of the FFLO phase. Previous proposals for the 
detection of the FFLO phase include detecting charac- 
teristic peaks in the expansion image following release of 



the trap [8|, the probing of noise correlations in an ex- 
pansion image [ij], and the use of RF spectroscopy to 
excite atoms out of the gas into other states [l^ • In a re- 
cent letter [l3] we showed how the FFLO phase in ID at 
T = can be detected via its collective mode signature. 
In this paper we extend this discussion to non-zero tem- 
perature and discuss other aspects of collective modes in 
imbalanced Fermi gases. Collective modes have proved 
very useful for probing t he p ro p erties of ultracold bal- 
anced [2I and imbalanced |l8l - [21| Fermi gases. Our work 
extends these discussions to one dimension and the pres- 
ence of the FFLO state. 

The paper is organised as follows. In Sec.|TT]we present 
the theoretical model of our study. In Sec. lIIII we describe 
the homogeneous system and present our results for the 
associated collective mode spectrum. We describe the 
speed of sound as a function of the polarisation of the 
Fermi gas and highlight the collective mode signature of 
the FFLO phase. We then discuss the effects of nonzero 
temperature. In Sec. IIVI wc turn our attention to the 
trapped Fermi gas. We assess the accuracy of the single 
mode approximation for various low-frequency modes in 
a trap. We then investigate the effect of the FFLO phase 
on the collective mode spectrum of the trapped imbal- 
anced Fermi gas. From this we deduce a signature of the 
presence of the FFLO phase in the collective mode spec- 
trum. Finally we describe the response of the trapped 
Fermi gas to a modulation of the coupling constant. 



II. THEORETICAL MODEL 

Wc study a model of a two-component Fermi gas in 
one dimension with attractive contact interactions and 
unequal densities. For a homogeneous system, the exact 
groundstates [1, [13, [22| - [23 | and thermodynamics [1^ are 
known from the Bethe ansat z l26l - i28| . As has been dis- 
cussed in detail by Liu et al. |29l|. Bogoliubov de Gennes 
(BdG) mean-field theory provides an accurate descrip- 
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tion of the exact phase diagram for the ID system, at 
least within the weak-coupling BCS regime. We have 
extended the mean-field theory to investigate the linear 
response and collective modes of both the trapped and 
untrapped imbalanced Fermi gas. 

Our numerical calculations are based on the discretised 
Bogoliubov de Gennes Hamiltonian which is given by 

+ E {^4Ai + ^•^•) + E W^,.cl,C.,a (1) 

where c-^j! are fermionic operators for species a =t, i on 



site «, W,^a = y^'^* + Vl 



Ha (with the external 
potential, V^^^ — U{c\gCi^g), /ig- the chemical potentials 
and a is the spin opposite to spin cr) and A.; = U {ci_^Ci_-f) 
is the local superfluid gap. J is the hopping parameter 
and U the on-site interaction strength {U < is assumed 
throughout). The angled brackets (...) denote the ther- 
mal and quantum expectation value at temperature T. 
The corresponding Bogoliubov de Gennes equations are 
solved self consistently to yield the mean field ground 
state. 

The results we present are at sufficiently low particle 
density to be representative of the continuum limit. In 
the mapping to the continuum, we relate site i to position 
X via X = ia, the mass is m = „" ;> and the ID interac- 
tion parameter is gid ~ Ua. The interaction strength 7, 
defined as the ratio of the interaction energy density to 
the kinetic energy density, is given by 



7 ^ 



t9p 



(2) 



[2 9( 1 where p is the total density of particles. In the case of 
a trapped system we take 7 to be defined by the density 
at the centre of the trap. We restrict our study to values 
of 7 up to about 7 = 1.6 where mean-field theory gives 
qualitatively accurate results [2^ . The Fermi wavevector 
for spin a is given hy kp„ = irp^ where p^ is the density 
for spin a. We define the polarisation p of the ID Fermi 
gas as 



(3) 



Additionally, we shall find it useful to define ksTp — 
Ep = TT^p^J/A. In the remainder of this paper we set 
h=l. 

The collective mode spectrum is obtained by consid- 
ering the linear response of a system to external per- 
turbations. Divergences in the response appear at the 
frequencies of the collective modes. We obtain the linear 
response of the equilibrium system by supplementing the 
Bogoliubov de Gennes equations with a self consistent 
random phase approximation (RPA). As was pointed out 



in and as we will argue in Sec. IIII CI this approxi- 
mate theory correctly captures all qualitative features of 
the collective modes that are important for our purposes. 



The self consistent linear response within RPA is com- 
puted as follows. Let 



(4) 



where pi,^ = c\,j{t)c^,„{t), k\{t) = c\.^{t)c\j^{t) in the 
Heisenberg representation. Furthermore, let Ai^a{i) (en G 
{1..4}) refer to the four components of Ai{t) and denote 
the Fourier transforms of Ai,a{t) as Ai,a{i^)- The bare 
response or single particle response of Ai^a to a pertur- 
bation SH{u!) = J2j p ^^],0i^)^j,l3i'^) is given by 



M,;,„(w)) =^xLj/3(w)5%^)(w) (5) 
where the bare response function x° is given by [soj 



xLj73(w) 



dt / dr 

— 00 J —00 



A,a{t),Ajp{T) 



Aujt 



(6) 



The thermal average (. . .) is computed by transform- 
ing to the quasiparticle eigenbasis in which (aj^a^) = 
Sfj.,iy.f{Efj,) where ajj is the quasiparticle creation opera- 
tor for a quasiparticle with energy E^ and / is the Fermi 
function. The self consistent response function x is ob- 
tained within RPA by taking into account the change in 

the Hamiltonian due to a change in the quantities (^Ai ui^. 

The perturbing potential will then consist of the external 
perturbation plus the potential arising from the interac- 
tion, 



(7) 



where |3C 



Hls'kj ^ U6i,kKs,^ and Ks, 
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Inserting this into the bare response we obtain: 

(m,,,(^)) = {k,5^.,-xl,isHirk,y\l,j0Shf; 

(9) 

(summation implied), which determines the RPA result 
for the susceptibility. 
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FIG. 1: (Colour online) Example of a homogeneous system 
displaying an FFLO phase. The periodicity of the densities 
is given by A a/2 where Aa is the wavelength of the order pa- 
rameter. This configuration has 7 = 1.6 and the polarisation 
p = 0.06. 



III. HOMOGENEOUS SYSTEM 

The ground state of a homogeneous Fermi gas with 
contact interactions can be found exactly by the Bethe 
ansatz [l^ . For attractive interactions the phase dia- 
gram spanned by the two chemical potentials at T = 
contains a fully paired region, a fully polarised (and 
therefore non-interacting) region and a partially po- 
larised region. The latter is associated with the FFLO 
phase. The mean field ground state was shown to be in 
goo d qualitative agreement with the exact ground state 
|29j . The FFLO phase is characterised by an oscillat- 
ing order parameter with a wavelength Aa and different 
average densities Po- for the t and l particles. As the 
pairing suppresses density differences, the density differ- 
ence p^{x) — Pi{x) is maximal at the nodes of the order 
parameter and minimal at the anti-nodes. For low po- 
larisations it is possible to describe the system as a fully 
paired supcrfluid on top of which there is a gas of excess 
majority spin particles. These are situated at the nodes 
of the order parameter [U [s^l ■ At zero temperature the 
wavevector = for the order parameter is given by 
kA = kp^-f — kpi- A possible ground state configuration 
where these quantities are illustrated is shown in Fig. [T] 
With increasing temperature the order parameter disap- 
pears. For interaction strengths we are considering, the 
FFLO phase disappears at temperatures T on the order 
oiT/Tp = 0.1 III. 



A. Sound velocities in an imbalanced gas 

We have investigated the speed of sound in a homo- 
geneous system and how it depends on polarisation. For 
this we compute the response of the system to potentials 




FIG. 2: (Colour online) Response of a homogeneous system in 
the FFLO phase to a spin-asymmetric periodic potential (|10[) . 
The area of the circle is proportional to the amplitude of the 
response. The polarisation is p = 0.15, and the interactions 
strength is 7 = 1.5. Two gapless sound modes are seen to 
emerge around A: ~ and fc ~ fc* = 2A:a. The simulation was 
done on 270 lattice sites. The divergent response at a; = is 
not shown. 



of the form 



5Wa {x,t) = Va sin kx cos u>t . 



(10) 



We refer to the case = Vi = Vo as a "spin-symmetric" 
and = — V4, = Vo as a "spin-asymmetric" excitation. 
Eq. pn]) represents a standing wave perturbation with 
a fixed phase. It is sufficient to consider this type of 
perturbation since we find that changing the phase of 
the perturbation with respect to the phase of the order 
parameter does not change our results in any significant 
way. Throughout this paper wc define the response R as 



R 



(11) 



which is a measure of the change in the density. 

For nonzero polarisations when the FFLO state is 
present two distinct sound modes exist with in general 
different velocities, as shown in Fig. [5] This can be un- 
derstood within mean field theory from the fact that the 
FFLO phase breaks both gauge and translational symme- 
try. However, as we will argue in Sec. IIII C[ these results 
are fully consistent with those expected from exact the- 
ory. The sound velocities at T = of these two modes as 
a function of the polarisation are shown in Fig. [31 (Small 
fluctuations in the lattice filling due to the techriique used 
lead to small variations in the sound velocity |34{.) For 
very low polarisations we effectively have a paired super- 
fluid on top of which there is a weakly interacting gas 
of excess majority spin particles, which are located at 
the pos itions of the domain walls in the order parameter 
[31I |32{. In this regime of low polarisation one of the 
sound modes is associated with phase fluctuations of the 
order parameter, and has a very similar velocity as that 
for a spin balanced superfluid. The other sound mode, 
with lower velocity, arises from the motions of the excess 
majority spin particles (the domain walls). The speed of 
sound of this mode becomes very small (below the reso- 
lution of our calculation) at very low polarisations, as the 
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FIG. 3: (Colour online) (a) Sound velocity versus polarisation 
for interaction strength 7 ~ 1.6. vp ~ ^{'VF,t + '''f,i)- As 
p ^ the velocity of one of the sound modes disappears. As 
soon as the domain walls between regions of positive A and 
negative A cease to overlap the sound velocity drops below our 
numerical resolution. This is because shifting the position of 
a domain wall costs very little energy, (b) Order parameter 
and density distributions for a polarisation (p = 0.038) at 
which the domain walls cease to overlap. 



domain walls become very dilute and cease to overlap. In 
the other limit, p — >■ 1, one of the sound modes vanishes, 
whereas the higher velocity mode tends to the value for 
a non-interacting Fermi gas. 



B. Signature of the FFLO phase 

We now turn to the signature of the FFLO phase in 
the collective mode spectrum. Again we investigate the 
response of the system to potentials of the form ([TUl) . 
However, this time we concentrate on short wavelength 
perturbations. 

For zero temperature there is a clear signature of the 
FFLO phase in the response spectrum of a homogeneous 
Fermi gas. As illustrated in Fig. [2] this manifests itself in 
sound modes emerging around k = k* = 2kA in addition 
to the ones emerging around fc = [iTj . Within mean- 
field theory, this can be understood as a consequence of 



the broken translational symmetry of the FFLO phase, 
though this too is fully consistent with exact theory (see 

Sec. Unci) . 

These low energy sound modes are due to a Brillouin 
zone (BZ) arising from the spatial periodicity of the su- 
perfluid order parameter A(x). The consequence is a 
generalised Bloch theorem. Let T be the operator which 
performs the transformation: x x -\- Aa/2, A(a;) — > 
— A(a;). For a system with an FFLO wavelength of Aa, 
T commutes with the Hamiltonian ([1]) since p-j- and pi 

have periodicities Aa/2. Note that T encodes the small- 
est translational shift for which such a symmetry exists. 
As the Hamiltonian commutes with T, the density matrix 
Pq for the equilibrium configuration is invariant under 
T^^PqT and T is time-independent. In the continuum 
limit the response of the system is given by 



(12) 

Within linear response is given by [35j 

X%ir,r',t,T) - i([i«(r,t),i;3(r',r)]). (13) 

As described above, in the self consistent random 
phase approximation SH is given by an externally 
applied perturbation plus a self consistent perturba- 
tion. Consider free oscillations, i.e. the external 
perturbation is set to zero. Then SH is given by 

dH{T) = I]/3 / dr' gidKp^^ (5A^{r\ r)^ with Kp^^ given 

by eq. ([5]). The equations of motion are then given by 



d_ 

di 



M„(r,t) 



+ J drj dr' ^^xlA^:r\t,T) 
guKp,J5A^{r\T)) (14) 



(summation implied). The equations of motion are in- 
variant under the operator T provided is invariant 
under T. This can be easily shown to be the case using 
the fact that T~^pqT = pq together with eq. p^ . 

The result of this generalised Bloch theorem is that 
fc* = (27r)/(AA/2) ~ 2k A is a reciprocal lattice vector 
for the collective mode spectrum. This means that the 
conserved quasimomentum is only defined up to multi- 
ples of fc*. Thus it is possible to perturb the system 
at a wavevector q and obtain a response at a wavevec- 
tor q — nk* {n integer) which allows low-frequency sound 
modes to be excited when q « nk* . Note that the gen- 
eralised Bloch Theorem only states that this is possible. 
The strength of this response needs to be calculated sepa- 
rately and is determined by the matrix element coupling 
perturbation q to response at q — nk*. A smaller am- 
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plitude of the oscillating densities, e.g. due to higher 
temperature, will lead to a weaker coupling of these re- 
sponses, see fig. m Once the system becomes homoge- 
neous the Brillouin zone structure disappears and the 
coupling vanishes. Note that this argument implies a BZ 
twice the size than if the periodicity was simply given by 
Aa. 



C. Mean field theory versus exact theory 

While the calculations we have presented are within 
RPA mean-field theory, as we now argue, the qualitative 
conclusions arc valid more generally. In the above, the 
response spectrum was accounted for in terms of the bro- 
ken translational and gauge symmetry of the mean-field 
state. As is well known, in a true ID quantum system 
no continuous symmetries are broken. Thus, the broken 
(phase and translational) symmetries of the mean-field 
FFLO state can lead only to power-law decay of the re- 
spective correlation functions [36j . However, the qual- 
itative features described above arc the same as those 
expected from an exact treatment of the system. The 
transition from the (unpolarised) superfiuid phase to the 
partially polarised phase is marked by the closing of the 
spin-gap, leading to a second gapless sound mode. This 
can be viewed as a Luttinger liquid representing the ex- 
cess fermions [3l|. These excess particles have density 
p-f — pi, and thus a Fermi wavevector {kp^^ — kF,i)- Thus, 
there are two gapless collective modes, arising from the 
fully paired particles, and the liquid of excess majority 
spin particles. Furthermore, as in the general theory of 
Luttinger liquids j37| the spectral function of these ex- 
cess majority spin particles will show gapless response 
at multiples of twice their Fermi wavevector. This is 
2(fc_F,-t' — kp^i) = 4:'K/Xi\ which is precisely the wavevector 
k* at which one finds the gapless response in mean field 
theory. Thus, the qualitative features of the collective ex- 
citation spectrum obtained in mean-field theory are fully 
consistent with those expected for the exact system. 



D. Temperature dependence of the collective 
modes 

We now consider the temperature dependence of the 
collective response of the homogeneous system. For very 
low temperature (T < O.OlTp) we find essentially the 
same response spectrum as for zero temperature. As 
the temperature of the system is increased and the order 
parameter disappears, the low energy response around 
k = k* also disappears, see Fig. U] This can be explained 
by the fact that the BZ associated with the oscillating 
order parameter disappears. We also observe that as 
the temperature is increased k* no longer coincides with 
2{kp^-^ — kF,i)- We note that thermally excited quasipar- 
ticles which do not contribute to the pairing can lead to 
this departure. 



The results we have presented so far are purely within 
the collisionless regime. The collisionless regime de- 
scribes the situation where the rate of quasiparticle col- 
lisions is much smaller than the frequency of the excita- 
tion. At zero temperature the system is always in the col- 
lisionless regime. On the other hand for non-zero temper- 
ature there are thermally excited quasiparticles. These 
may undergo collisions, thereby leading to a nonzero col- 
lision rate. As pointed out in Ref. [33], at nonzero tem- 
perature great care has to be taken when taking the limits 
fc — >■ 0, T — >■ 0. These two limits do not commute. If first 
the temperature is taken to zero and subsequently k is 
taken to zero, the system is in the collisionless regime. 
Sound modes in the collisionless regime correspond to 
zero sound. On the other hand if first fc — > and sub- 
sequently T — > a hydrodynamic description is needed. 
In the hydrodynamic regime quasiparticles have time to 
scatter many times during the cycle of one collective oscil- 
lation and hence the system can be described in terms of 
local thermodynamic equilibria. In the crossover region 
no sound propagation is possible due to strong damping. 

We have performed a simple analysis of the quasipar- 
ticle collision rate. In an infinite one-dimensional system 
no quasiparticle excitations can exist due to the diver- 
gent nature of the interactions; a Luttinger liquid [37| 
description of the system is required [39j . All excitations 
are collective excitations. This manifests itself in an infi- 
nite quasiparticle scattering rate. On the other hand, in 
a finite system, the quasiparticle scattering rate is finite. 

Using perturbation theory we find that in a system of 
size L, for T <^ Tp, the scattering rate of a quasiparticle 
with energy ksT is 



SnhEp e + 1 



T 



2 In 



An Tp' 



kpL 



(15) 



In general for an infinite system the collisionless ap- 
proximation will fail since F diverges. In this case a hy- 
drodynamic theory must be used. However, in practice 
all experimentally relevant systems are finite. There- 
fore there will be some temperature below which the 
system can be described as collisionless. For example, 
comparing the frequencies with the quasiparticle scat- 
tering rate for the system sizes used to compute Fig. SI 
we find that in all apart from the four lowest frequency 
points for T = O.OTGTp and the lowest frequency point for 
T = 0.057Tp the systems are in the collisionless regime. 
We therefore expect the results presented in Fig. 0] to be 
qualitatively accurate for systems of this size. 

For a given system of size L the temperature below 
which the system can be considered collisionless may be 
found as follows. The lowest k vector with which the 
system can be excited is given by fc = Fig. [3] suggests 
that for the parameters we are considering, the sound 
velocity v is of the order of the Fermi velocity vp. Then 
using the fact that the collisionless approximation will 
be valid up to approximately w w F and using uj « kvp 
we obtain that for a homogeneous system the maximum 



6 




^^1 Iti T^io ^Sm It5 T^iD oTi Jl3 T^S Tjs 150 HT^ ^tS u 

fc/(2fci) k/i2k^) k/(2k^) k/i2k^) 



FIG. 4: (Colour online) Response of a system at fixed chemical potential and different temperatures. The area of the circle is 
proportional to the amplitude of the response, as in Fig. [J] The polarisation is p = 0.15, the interactions strength is 7 = 1.5 
and there are roughly 135 particles. The temperatures are from left to right T/Tf = 0.0037,0.038,0.057,0.076. Within our 
mean field theory calculation the oscillating order parameter disappears at a temperature T /Tf = 0.07. Simultaneous with the 
disappearing of the oscillating order parameter the response around k = k* disappears. 



temperature for which the coUisionless approximation 
is vahd is given by 

Here we we have only included the logarithmic depen- 
dence of r on T /Tp and not the linear part since the 
linear part is much smaller. Using cq. ([T5|) we find that 
for the configurations in Fig. 2] the coUisionless approxi- 
mation is valid up to temperatures of around T = O.OSTi?. 
For these systems this temperature coincides roughly 
with the temperature for which the FFLO phase disap- 
pears. 

Above we argued that for the systems and tempera- 
tures studied in this paper we can work in the collision- 
less regime. The results of our calculations are valid in 
the regime in which interactions are assumed to be weak 
or moderately strong. On the other hand it has been re- 
ported that for strongly interacting 3D imbalanced Fermi 
systems one needs to consider the effects of collisions 
of quasiparticles even at experimentally achievable low 
temperatures [l^, . There are some fundamental dif- 
ferences between those studies of 3D systems and the 
ID systems of interest here. In the 3D systems, a spin 
imbalanced phase is stable only at relatively high po- 
larisations (beyond the Chandrasekhar-Clogston limit). 
The collision rate is rapid until the temperature is small 
compared to the Fermi temperature of the minority com- 
ponent (40j . which, in view of the large imbalance, can 
be a low energy scale. In contrast, in ID the FFLO 
phase occupies a large region of parameter space, and it 
is possible to have this spin-imbalanced phase with much 
smaller polarisations when the Fermi energies for spin 'f 
and \, atoms are of the same magnitude. One then ex- 
pects that, for strong interactions, the relevant energy 
scale below which collisional scattering starts to become 
small is just the (typical) Fermi temperature. Our formu- 
lae (fTS]) and (fTH]) were derived using perturbation theory, 
so cannot be applied to determine the collisional scatter- 



ing rate in the strongly interacting regime. We are not 
aware of a calculation of the thermal broadening of the 
collective modes for a strongly interacting spin-i Fermi 
system of relevance here. 



IV. TRAPPED SYSTEM 

We now turn our attention to the experimentally 
more relevant case of the Fermi gas in a harmonic trap. 
Within the local density approximation (LDA) the pos- 
sible ground state configurations for ID Fermi gas at 
T = are known exactly from the Bethe ansatz [l^ . 
There are only two possible configurations which are (a) 
a partially polarised phase in the middle, with a fully 
paired phase towards the edge and (b) a partially po- 
larised phase in the middle, with a fully polarised phase 
towards the edge. Recently these configurations have also 
been observed in experiment [l^ . As shown by Ref. [2^ 
Bogoliubov de Gennes theory reproduces these ground 
state configurations and it is in good qualitative agree- 
ment with the results of the exact solution combined with 
the LDA. 

The low-frequency collective modes in a trap are the 
dipole and the breathing modes. Dipole modes are the 
lowest lying modes with odd parity around the centre 
of the trap. For T4xt(a;) = ^mwgX^ their frequencies in 
the non-interacting limit are given by w = cjq. Breath- 
ing modes are the lowest lying modes with even parity. 
They correspond to expansion and contraction of the 
atomic cloud. In the non-interacting limit their frequen- 
cies are given by a; = 2aJo- Both the dipole and breathing 
modes exist in two types: the spin and density modes. In 
the density modes the two atomic clouds move in phase 
whereas in the spin modes the two clouds move largely in 
anti-phase. Within these modes the density dipole mode 
has a special status. It is called the "Kohn mode" and 
corresponds to rigid body oscillations of the cloud about 
the centre of the trap. Its frequency uj is given by w = 
independent of the interaction strength. 
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We now investigate the accuracy of the single mode 
approximation for these low frequency modes. 



A. Single mode approximation 

The single mode approximation (SMA) is an important 
tool in calculating collective mode frequencies. Together 
with sum rules it allows the calculation of collective mode 
frequencies from ground state properties alone. This can 
be done as follows: for an operator F the dynamic form 
factor S{ui) describes the strength of the transition from 
the ground state to state with energy huj. S{uj) is given 

by m 



5M -^|(fc|i^|0) S{uj~LOk) (17) 

where |fc) label the exact energy eigenstates of the sys- 
tem. We define the pth moment of the dynamic form 
factor as 



tp — I S{(jj)uj^ d(jj . 
/o 



(18) 



Using sum rules it is possible to express mi, and m_i 
as Ell 



nil 
ma 

m_i 



H.F 



F,H 



H,F 



-lmx{F,LO^ 0) 



(19) 
(20) 
(21) 



where xi^^^) is the response function for the operator 
F at frequency tu. 

The SMA corresponds to the assumption that S{u!) is 
dominated by one sharp peak. Then the collective mode 
excited by the operator F has a frequency that is given 

hy Up ^ i^p,p-2 = \l m^'' 2 ^'^"^ value of p. If the SMA 
is not exact, and more than one peak is present, one 
finds that Wp.p_2 > UJq^q-2 for p > q and this approach 
to finding the collective mode frequency using sum rules 
becomes less accurate. Nevertheless, a;p,p_2 still provides 
an upper bound for the lowest frequency collective mode. 
We use both and uj^^i to obtain upper bounds for 

the collective mode frequencies. 

In order to investigate the dipole and the breathing 
modes we consider two sets of operators F. For the 
dipole mode we consider Fd{0) = cos 6 J2i a^i +sin 9 J^j Xj 
where indices i run over all the t particles and j over all 
the 4, particles. For the breathing modes we consider 



Fbie) 



sin J2j ^'j ■ Pd corresponds to ap- 



plying hnear potential V^{x) = xcosO ,Vi{x) = xshid 
and Ff, corresponds to applying a quadratic potential 



V^{x) = x^ 



• Viix) 



Here the single mode approximation amounts to as- 
suming that the operator F{9) has at most two peaks 
in the response spectrum, each of which can be made 
to vanish by suitable choice of 6. Hence the collective 
mode frequencies of the two associated modes can be ob- 
tained by maximising or minimising ujp_p^2- The lower 
frequency mode for each F is the density mode and the 
higher frequency mode is the spin mode. [48| 

For Fd we obtain in the continuum limit two values 
of 9 which maximise or minimise CJ3 1. For both T = 
and for T > these are 9sd = — arctan for the spin 
dipole mode and 6'Kohn — f for the Kohn mode. For the 
breathing mode such simple analytical formulae for 9sb 
(spin breathing mode) and 9^,1, (density breathing mode) 
which maximise/minimise ujp^p-2 could not be obtained. 
These can be obtained numerically though as has been 
done in Ref. [lo| . where the density breathing mode fre- 
quency was investigated using sum rules. 

We have used our results to address the question: how 
accurate is the SMA in describing the results of the 
BdG+RPA approximation? 

For all the modes mentioned above we find that the 
SMA is strictly speaking not correct, since uj^^i 7^ 
To gauge the accuracy of the SMA we have compared the 
prediction of the sum rule for with the position of 

the peak in the response spectrum, see Fig. [5l Wc choose 
uji^-i because it gives higher accuracy results as it is not 
affected as much by low weight high frequency peaks as 
Lj^^i is. For the two density modes which are the lower 
frequency modes of the two operators F^ and Ff, we find 
that there is very good agreement of the positions of the 
peak as obtained from the response spectrum and the 
SMA. For the density dipole mode (the Kohn mode) the 
sum rule calculation captures the position of the peak 
up to our numerical accuracy. For the density breathing 
mode the accuracy typically is about 1%. In the limit 
p ^ 1 the single mode approximation becomes exact for 
these two modes as can easily be verified by evaluating 
eqs. (fT9|) and ((20|) for a non-interacting gas. 

However, for the spin modes the SMA is not as accu- 
rate. In particular for lower polarisations the frequency 
from the SMA estimate coi^-i can differ quite significantly 
from the position of the lowest frequency peak of the full 
spectrum. This is due to the presence of many small am- 
plitude high frequency peaks. The increased weight of 
these high frequency peaks and the associated failure of 
the SMA is correlated with an increase in size of the fully 
paired region towards the edge of the cloud. 



B. Collective mode signature of the FFLO phase 

As described in [l3| at T = we find a dramatic sig- 
nature of the FFLO phase of a trapped Fermi gas when 
the spin-dipole mode is excited by a short-wavelength 
potential. When perturbing with potentials of the form 
pop we find a large response of the spin dipole mode 
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FIG. 5: (Colour online) Calculation of the accuracy of the 
SMA using (a): breathing modes and (b): spin modes. 

These plots are for 7 — 1.2 and different polarisations. For 
the density modes (a) the single mode approximation is very 
accurate. For the spin modes (b) the single mode approxi- 
mation becomes inaccurate at small polarisations. The arrow 
indicates the point where the outside of the cloud goes from 
being fully paired (p < 0.1) to being fully polarised {p > 0.1) 
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FIG. 6: (Colour online.) Configuration with fully paired 
phase towards the edge of the trap, (a); densities and the 
value of the superfluid gap A. (b): response of the spin- 
dipole mode to excitations of different wavelengths (arbitrary 
units). Here p — 0.048, 7 — 0.93 (measured in centre), num- 

_o 1 

ber of particles A*' — 290, lattice spacing a — 3.3 • 10 N2 auo 
with aho = •\J^^~- The perturbing potential has a fixed am- 
plitude, while the wavelength is varied. 



when the wavelength of the perturbing potential becomes 
Aa/2 [131 . This response is independent of whether 
the ground state configuration is of type (a) or (b) , that 
is whether the outside is fully paired or fully polarised. 
This response for these two types of system is shown in 
Figs. [S] and [71 This response can be understood in 
terms of coupling to the gapless modes at A: = 4^- = fc* 




FIG. 7: (Colour online.) Configuration with fully polarised 
phase towards the edge. Here p = 0.25, 7 = 1.5, a/ {N^aho) = 
4.7 ■ 10"^, = 143 . Otherwise the same as Fig.[ni 



of the homogeneous system (see Figl2]and Sec. IIII B|) . A 
small deviation from k* (of the order the inverse system 
size) allows mixing of this mode to the spin-dipole mode, 
causing the response at k* to be apparent in the dipolar 
motion of the atomic cloud. 

As in the untrapped system, at very low temperatures 
the system behaves in the same way as for T = 0. With 
increasing temperature the oscillating gap A disappears. 
Along with it the sharp response at the wavelength Aa/2 
disappears, see Fig. [H There are a few additional note- 
worthy points. As the temperature is increased the num- 
ber of modes which can be excited increases. This can 
lead to the splitting of, for example, the spin dipole mode 
into two near degenerate sub modes. For yet higher tem- 
perature further fractionalisation of the modes occurs, 
see Fig. [Sj At intermediate temperatures (Fig. 
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oscillating gap can still be made out but the peak in the 
response of the spin dipole mode is strongly attenuated, 
which may increase the difficulty of finding the FFLO 
phase. An example of a system which at T = is fully 
paired towards the outside is given in Fig. HI but the case 
where the outside is fully polarised at T = is qualita- 
tively similar. 

As in the homogeneous case discussed in Sec. IIII B] our 
calculation is performed in the coUisionless regime. Us- 
ing eq. ((T6)) and the LDA by taking the system size to 
be roughly the size of the cloud we find that the coUi- 
sionless regime should be valid up to a temperature of 
T/Tf « 0.03. This is slightly lower than in Sec. HITdI 
since the system size is larger. Nevertheless, the three 
configurations in Fig. [5] with an FFLO phase are still in 
the coUisionless regime. 

Most results presented in this paper are for values of 
7 w 1 — 1.6 where the mean- field theory we use is expected 
to be accurate. (For the groundstate, this is confirmed by 
comparisons with exact results [1^.) We have performed 
calculations with larger 7, in the range of 7 = 2.2 — 3.8, 
and find qualitatively the same results. While we can- 
not trust the quantitative results of mean-field theory in 
regimes of strong interactions, particularly at unitarity, 
we believe that the qualitative results remain valid. In 
particular, the key signature of the FFLO phase that we 
predict relies only on the commensurability between the 
wavelength of the perturbing potential and the intrinsic 
periodicity of the FFLO phase, and is expected to be 
robust. 



C. Experimental considerations 

Recently experiments have produced true ID ultracold 
Fermi gases as an array of ID tubes [l^. By changing 
the coupling between the tubes it is possible to tune con- 
tinuously between a 3D system and a true ID system. 
The experimental results show evidence for the appear- 
ance of a partially polarised phase in a regime of pa- 
rameters consistent with the theoretical expectations for 
the phase associated with the FFLO state. Probing the 
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FIG. 8: (Colour online) Response of a system (for which at T = the outside is fully paired) to a periodic potential of 
wavelength A and freqency ui. There are 173 particles, 7 = 0.93, p varies from 0.05 to 0.09. The temperature is T = O.OOOSTf, 
T — 0.005Ti?, T = 0.02Tp and T = .IOGTf. For the lowest temperature the spin-dipole frequency is oi = 1.31ajo. The sharp 
peak at the spin-dipole frequency at \/ {N^^^ aho) ~ 0.2 is the signature of the FFLO phase. At higher temperatures, this sharp 
response disappears and additional excitations emerge. The case where the outside is fully polarised towards the edge at T = 
is qualitatively the same. 



collective modes in the manner described in this paper 
would lead to direct evidence of the intrinsic modulation 
of the FFLO phase. 

Observing the signature of the FFLO phase in experi- 
ments in this way requires the ability to create a variable 
wavelength optical lattice. This can be done using the 
technique used by Ref. [l^. Although the signal ap- 
pears for a spin-symmetric perturbation, the signal is 
clearer for spin-asymmetric perturbations. Thus, any 
spin-dependence of the optical lattice will enhance the 
ability to distinguish the peaks associated with the os- 
cillation of A from the background peaks. For Bosons 
(^^Rb) potentials with spin-dependence have been cre- 
ated by Ref. jU, and would be ideally suited to this 
purpose. There are additional complications associated 
with applying spin-dependent potentials to ^Li due to 
increased loss rates resulting from the application of the 
potentials [i^l ■ As discussed in Ref. Q it is nevertheless 
in principle possible to set up these potentials. There are 
two natural ways to find the response of the spin-dipole 
mode. One way is to follow the approach of the calcula- 
tion precisely, and make the strength of the optical lattice 
time dependent, as in Ref. [i^. The spin-dipole mode 
can then be selectively excited by bringing the temporal 
oscillation of the lattice in resonance with the spin-dipole 
mode. Alternatively a static periodic potential can be 
applied to the system and the system allowed to equili- 
brate. Subsequently this potential is abruptly switched 
off (similar to the approach used in Ref. ^^). This will 
excite collective modes of many frequencies from which 
the response of the spin-dipole mode can be obtained by 
Fourier transform. 

In it is maintained that in order to see the FFLO 
phase by in situ imaging it is advantageous to induce 
coupling between the different ID tubes in order to align 
the phases of the FFLO state between the different tubes. 



However, for the method we are considering this should 
not be necessary. The response of the system depends on 
the wavelength with which the system is excited and on 
Aa- Only one laser is used to create the optical lattice 
with which the system is excited which implies that the 
wavelength with which the system is excited is the same 
for all tubes. Therefore, as long as the wavelength of the 
FFLO phase is approximately the same between different 
tubes, the response will be be small or large for the same 
excitation wavelengths. As a result it should be possible 
just to work with very weakly coupled or uncoupled ID 
tubes and still see this effect. 

One experimental difficulty in detecting the FFLO 
phase in systems of the type studied in Ref. - which 
is shared by all probes of the intrinsic periodicity of the 
FFLO phase - is the inhomogeneity across the many ID 
tubes. Due to the different positions of the tubes in the 
atomic trap, different tubes will have different chemical 
potentials. In general this will lead to different values 
of Aa in different tubes. Since the large response occurs 
for perturbations with A = Aa/2 only those tubes which 
meet the condition on Aa will show a strong response. If 
all the tubes can only be imaged as a whole it is therefore 
desirable to work in a regime where Aa varies minimally 
between different tubes. As pointed out by Ref. [l3 it 
is possible to choose a point in the phase diagram where 
= = 0. At this point the value of A a should 

be the same over a large range of tubes in the centre of 
the system. 

Our method provides a very elegant way to overcome 
the effects of inhomogeneity, if in situ imaging can allow 
individual or a small collection of ID tubes to be resolved. 
This will allow one to relax the condition of having any 
intertube tunneling and being at a special point in the 
phase diagram. By exciting the whole system with a 
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given wavelength A, only those tubes with Aa — 2A will 
have their spin-dipole oscillations excited. Tubes with 
the same value for Aa lie on concentric cylinders around 
the trap centre. Thus, if imaging of the cloud has suf- 
ficient resolution to identify which tubes are oscillating, 
one should see response from tubes that lie on a cylin- 
drical region around the axis of the trap. Changing the 
wavelength A will change the radius of this cylindrical re- 
gion. In this way the FFLO wavelength could be studied 
as a function of the position of the ID tube in the trap. 
This approach is only limited by the power to resolve the 
ID tubes. 



D. Modulation of the coupling constant 

We have investigated the response of the one- 
dimensional Fermi gas to a modulation of the coupling 
constant. This is an interesting perturbation to consider, 
since using Feshbach resonances it is relatively easy to im- 
plement a uniform modulation of the coupling constant 
experimentally. Recently it was also found, that in a lat- 
tice system parametric modulation can be achieved by 
modulating the lattice depth j47| . 

Modulations of the coupling constant couple to den- 
sity differences. The FFLO state has a density differ- 
ence that shows quasi-periodic spatial oscillations, as ex- 
plained above. Hence it is interesting to ask if the re- 
sponse of the system is sensitive to the appearance of the 
ordering characteristic of the FFLO phase. 

Within linear response the modulation of the coupling 
constant of the form U = Uo + SU (t) leads to a pertur- 
bation of the form 



(22) 
(23) 



and hence the response to a modulation of the coupling 
constant is amenable to our RPA calculation. 

We have studied the response of the trapped gas to 
modulations of the coupling constant, focusing on the 
regime of low frequencies. Since a modulation of the cou- 
pling constant does not break refiection symmetry about 
the centre of the trap it can only excite trap modes of 
even parity. As a result the lowest frequency modes which 
can be excited arc the breathing modes. 

We find that the nature of the response of the system 
depends strongly on whether the system is in the con- 



figuration (a), which is fully paired toward the edge and 
exhibits an FFLO phase in the centre, or in the config- 
uration (b), which is fully polarised towards the edge. 
When the outside is fully polarised the majority of the 
response weight is at the lower frequency density breath- 
ing mode and there is little or no response at the higher 
frequency spin breathing mode. On the other hand, in 
the case where the outside is fully paired and there is an 
FFLO phase in the centre the situation is reversed. Most 
of the weight of the responses on the high-frequency spin 
breathing mode. See Fig. O 

While the nature of the response is senstive to the 
global distribution of particles, whether or not an FFLO 
phase is present makes little difference to the response 
spectrum. We must therefore conclude that there is no 
clear signal of the presence of the FFLO phase in the 
collective mode spectrum excited by a modulation of the 
coupling constant. On the other hand the collective mode 
spectrum appears to be sensitive to whether the outside 
region of the cloud is fully paired or fully polarised. 



V. CONCLUSION 

To conclude, we have investigated the collective mode 
spectrum of the imbalanced Fermi gas. The FFLO phase 
can be detected due to the formation of a characteristic 
periodicity of the densities. This makes it possible to use 
short wavelength perturbations which couple to the in- 
trinsic periodicity of the FFLO phase in order to excite 
long wavelength responses which are easy to measure in 
experiment. This signal in the collective mode spectrum 
persists up to the temperature where the FFLO phase 
disappears. We also investigated the accuracy of the sin- 
gle mode approximation. We found that it works very 
well for density modes, but not as well for spin modes. 
Finally, we found that a modulation of the coupling con- 
stant results in a collective mode spectrum which allows 
us to distinguish long wavelength characteristics of the 
density profile but not short wavelength characteristics, 
as are associated with the FFLO phase. 
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